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ABSTRACT 

The elastic anomaly observed in the coherent Kondo state of Ce heavy 
fermion compounds is analyzed by using the Anderson lattice model sim- 
ulating the energy level scheme of CeTe. The T 7 doublets and T 8 quartets of 
the 4f states are considered in the model. We solve the mean field equations 
to derive the temperature dependences of elastic constants, using the random 
phase approximation like expression for the interaction between the elastic 
strain and the crystalline field splitting. We compare the calculation with 
the (cxi — ci 2 )/2 and c 44 constants of CeTe. The presence of the downward 
dip and the observed overall temperature variations of the two constants are 



well described by the present theory. The origin of the dip is the coupling 
between the elastic strain and the splitting of the T 8 quartets. 
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I. INTRODUCTION 



Heavy fermion compounds show elastic anomalies in low temperatures. They 
are related with crystalline field structures, magnetic phase transitions, and 
so on. The main contribution coming from the crystalline fields is the peak 
or dip structures in the temperature dependences of the elastic constants. 
Most of the crystalline field splittings are larger than the Kondo temperature 
and the anomalies occur in higher temperatures than those of the coherent 
Kondo state. For example, the elastic constant c 33 of CeCu 6 [1] has a dip at 
about 10K. This is due to the splitting larger than the Kondo temperature 
4K. However, a few compounds have the splittings which are comparable 
to or smaller than the Kondo temperature. The anomalies occur in the 
coherent Kondo state. The constants (en — ci 2 )/2 and c 44 of CeTe [2] show 
the apparent dip at about 15K. The j = 5/2 levels of Ce ions split into the 
T 7 Kramers doublet and the T 8 quartet states. The T 7 states are the ground 
states. There is the splitting 30K between T 7 and T 8 states. This is the origin 
of the dip. Similarly, the c 44 constant of the alloy Ceo.5Lao.5Bg shows the dip 
at about 0.2K, and this temperature is lower than the Kondo temperature 
of this alloy [3]. 

The main purpose of the present paper is to develop a theoretical de- 
scription of the elastic anomaly by using a microscopic model. We use the 
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Anderson lattice model simulating the energy level scheme of CeTe. The I\ 
doublets and T 8 quartets are considered. We solve the mean field equations 
and assume that a strain field couples linearly with the splitting of the T 8 
states. We calculate the linear susceptibility with respect to the strain field. 
Then, we derive the temperature dependence of the elastic constant, using 
the random phase approximation (RPA) like expression [4] which includes 
the elastic susceptibility. We compare the calculation with the (en — ci 2 )/2 
and C44 constants of CeTe. The observed overall temperature variations of 
the two constants are well described by the present theory including the 
coupling between the elastic strain and the splitting of the T 8 quartets. 

We explain the model in Section II. We report the solution in Section III 
and the elastic properties in Section IV. We summarize the paper and give 
discussion in Section V. 



II. FORMALISM 

We formalize the infinite-f/ Anderson lattice model simulating the crystalline 
field structures of CeTe [2]. We use the slave boson method. The model has 
the following form: 

H = £[#f £+(£f + A - $ £+(£ f + A + 5) Y,\flki (!) 

i 1=1,2 1=3,4 1=5,6 
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+ E £ kC^C ki , 
k,/=l-6 

i,/=l— 6 

i /=l-6 

where is an annihilation operator of the f-electron of the l-th orbital at the 
i-th site, Ck,« is an operator of the conduction electron with the wave number 
k, and bi is an operator of the slave boson which indicates the unoccupied 
state at the f-orbital. The atomic energy of the first and second orbitals of 
f-electrons is Et\ that of the third and fourth orbitals is E{ + A — 5; and 
that of the third and fourth orbitals is E { + A + 5. The two crystalline field 
splitting parameters A and 5 are considered in the model. The first one A is 
the splitting between the ground state T 7 doublet and the excited T 8 quartet 
states of f-electrons. The second splitting 5 is due to the lattice distortion 
from the cubic symmetry. It is assumed that 5 couples with a strain field 
e linearly: 5 = r/e, where i] is the coupling constant. For the conduction 
electrons, the same quantum number is assumed as that of the f-electrons. 
We use the square density of states, p = 1/ND, which extends over the 
energy region, — D < < (N — 1)D, where N = 6 is the total number of 
quantum states. This assumes that the combination NpV 2 , which appears 
in the 1/N expansion, is independent of N. Therefore, the mean field theory 
becomes exact as N — > oo. The third term in the hamiltonian is the mixing 
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interaction between f- and c-electrons, V being the interaction strength. The 
last term limits the maximum number of f-electrons per site up to unity. 
This could be realized by the constraint J2i=i-6 n,ifi,i + b\bi = 1 with the 
Langrange multiplier field Aj. 

This model is treated within the mean field approximation: = r, 
(b\bi) = r 2 , and Aj = A (a site independent real value). These mean field 
parameters are determined by solving the following coupled equations [4]: 
(1) the constraint condition, 



3D J 



+ — / dE^ f(E - n) 

3D J (E { + A-5-E) 2 V W 

+ — [ dE^— f(E-n) +r 2 = 1, 

3D J (£f — E) 2 ' 

(2) the self-consistency condition for r, 



2 



-/ 

3D J 

lr V 2 

+ 37j/ d£ ^ /(£ -" ) + A = ' 

and (3) the conservation condition of electron number n c \, 

wI^ 1+ { E, + ^S-Ef WE -^ (4) 
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where f(x) = l/[exp(x/T) + l] is the Fermi distribution function, E { = Ef+\ 
is the effective f-level, and V = rV is the effective mixing interaction. The 
integrations are performed over all the energy region of the bands. The three 
equations are solved numerically for the three variables, r, A, and the Fermi 
level fi. In addition, the values at T = can be obtained analytically. 



III. SOLUTION 

Equations (2), (3), and (4) are solved numerically for the parameters D = 
5 x 10 4 K, V = 7500K, E t = -10 4 K, and n cl = 1.9 as the typical values. 
We take the splitting parameter A = 30K. The close value for A has been 
revealed in the experiment [5]. We consider the limit 5-^0 because the 
splitting coming from the strain is so small and negligible. As we will see 
later, A is smaller than the Kondo temperature T K = Ef — fi at T = OK. 

Figure 1 shows the temperature dependences of parameters. Figures 1 
(a), (b), and (c) show the variations of Ef, T K , and the number of f-electrons 
per site rif, respectively. As the temperature increases, the order parameter 
r decreases, so that n f = 1 — r 2 increases. The quantity r does not vanish 
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even though the temperature is much higher than 7k (about 40K) at T = 0. 
This is the effect of the change of the Fermi level /x to keep the total electron 
number constant. This effect has been reported previously [6,7]. According 
to the increase of rif, Ef decreases, which means the reduced itinerancy of 
f-electrons owing to the increase of rif. At low temperatures, the excitation 
energy is limited by the smaller distance from the Fermi level to the gap 
of the bands I = 1,2. This results in the increased value of rif when the 
crystalline field is switched on. Also, E { decreases and T K increases, due to 
the crystalline field. The similar dependence on A has been reported in the 
previous paper [4]. The Kondo temperature 7k (A) (at T = 0) as a function 
of A satisfies the equation, [7k (A) + A] 2 T|(A) = T|(0), where T K (0) = 
Dexp[— D((j, — Ef)/V 2 } is the Kondo temperature for A = 0. Starting from 
this analytic expression, we could verify the low temperature variations of 
parameters by using the expansion with respect to A/T K (0) assuming the 
small A. 
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IV. ELASTIC ANOMALY IN LOW TEM- 



PERATURES 



We shall discuss the change of elastic properties of heavy fermions due to the 
crystalline field splitting in the low temperature below T K . We shall calculate 
an elastic constant c by the RPA-like formula [4] analogous to the plasmon 
excitation theory. The constant c is related with the linear susceptibility 
with respect to 5, as shown below: 



where Co is the elastic constant of the system where there is not interactions 
between the lattice and the electronic system, and g is the coupling constant. 
There is a relation g = coi] 2 , so g is positive. The analogous formula was 
used before [8] but in the linear response theory. We assume that c is 
independent of the temperature. The value of g is unknown experimentally 
as well as theoretically. In order to discuss the crystalline field effect on c, 
we treat the factor g as a kind of fitting parameters. The quantity xs is 
calculated as the second order derivative of the mean field free energy: 



Xs = 



d 2 F 



(6) 




V 2 



(Ef + A-5-E) 
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V 2 



+ 7777 



(£f + A + 5 - E) 



where the S dependences of the band edges are neglected in the derivatives 
because their effect is exponentially small. In the actual calculation, we 
take the limit 5^0, because our problem is the elastic property at the 
equilibrium of the cubic lattice where there is not the splitting S. 

Figure 2 displays the temperature dependence of Xo = ^ m s^oXs- Figure 
2(a) shows the variation over wide temperatures, and Fig. 2(b) shows the de- 
tailed structure in low temperatures. There is a peak around T = 15K. This 
is owing to the large degree of freedom for electrons and the crystalline field 
A. The position of the peak would depend on parameters, but here the posi- 
tion agrees with that of the (en — ci 2 )/2 constant of CeTe [2]. The appearance 
of the peak has been discussed in the previous paper [4]. The susceptibility 
in high temperatures is nearly inversely proportional to T, showing the Pauli 
paramagnetic behavior. The value of Xo at T = is analytically expressed 
as, 



by using the Kondo temperature. 

Now, we compare the calculated c/co with the experiments. We plot 
the temperature dependences c/c , which are obtained from the experimen- 
tal data of CeTe [2]. They are shown by the dots. Figures 3 (a) and (b) 



2T K (A) 
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Xo = 



[T K (A) + A][3T K (A) + A]' 



are for the (en — ci 2 )/2 and C44 constants. The experimental Co depends 
on the temperature. We used the linear dependence used in Ref. 2 for the 
(cu — ci2)/2 mode. The quantity Co becomes softer as the temperature in- 
creases. However, we cannot use the experimental Co for the C44 mode [2] in 
order to compare with the theory. We rather use the increasing linear func- 
tion: c (T) = 0.696 + 9.89 x 10- 5 T(10 n erg/cm 3 ). It seems strange that the 
constant becomes larger as T rises. But, this does not mean that the crys- 
tal becomes harder for increasing T. Several lattice constants can become 
harder as T increases. In fact, c 44 becomes harder, and Cn and (en — Ci 2 )/2 
become softer, in CeTe. In the two figures, the calculated c/co is shown by 
the curve for g = 6.9K. The elastic constant decreases from much higher to 
lower temperatures than Tk. The decrease is almost proportional to 1/T. 
The agreement is good enough. The decrease is the effect of the valence fluc- 
tuation. There is a downward dip around 15-20 K and the position agrees 
with the experiments, too. The overall temperature dependences are well 
explained by the same g for the two constants. Of course, both Co and r] 
are different for the two constants. Thus, this result should be regarded as a 
coincidence. 
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V. SUMMARY AND DISCUSSION 



We have solved the mean field equations of the Anderson lattice model with 
the crystalline field splitting between T 7 doublets and T 8 quartets. It has 
been assumed that the strain field couples linearly with the splitting of the 
T 8 states. We have calculated the linear susceptibility with respect to the 
strain field. Next, we have derived the temperature dependence of the elastic 
constant, using the RPA-like expression which includes the elastic suscepti- 
bility. We have compared the calculation with the (en — Ci 2 )/2 and c 44 
constants of CeTe. The observed overall temperature variations of the two 
constants are well described by the present theory including the coupling 
between the elastic strain and the splitting of the T 8 quartets. We believe 
that the presence of the peak in xo is n °t an artifact of the mean field theory. 
In fact, the magnetic susceptibility of the exact solution of the single site 
system has a peak when the number of the degree of the freedom is larger 
than two. 

In the actual compound, V is anisotropic, i.e., it has an angle dependence 
in the momentum space: V = V{6,4>). The mean field equations, Eqs. (2), 
(3), and (4), change only in one point: the angle average, 



appears in the equations. This angle average can be absorbed in the present 
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formalism by regarding the angle averaged as the isotropic V of Eq. 

(1). Thus, the parameter value V used in this paper should be interpreted 
as an averaged one. 

In Ref. 3, the alloy system Ceo.5Lao.5Bg has the downward dip in the 
temperature dependence of the C44 constant. The variation is very similar 
to that in Fig. 3. The same mechanism of the elastic anomaly discussed 
in this paper would work in this alloy system, too. The elastic anomaly 
in magnetic alloys could be treated by a microscopic theory by using the 
coherent potential approximation applied to the Anderson alloy system [9]. 
This calculation will be an interesting extension of this paper. 

By applying the magnetic field, the atomic energy levels of Ce will split 
further. This will result in more structures in temperature variations of 
the elastic constants. Experimental as well as theoretical information will be 
useful for detailed understanding of the electronic properties of heavy fermion 
systems. 
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FIGURE CAPTIONS 



FIG. 1. Temperature dependences of the mean field solution: (a) E { , (b) 
T K , and (c) n { . Parameters are D = 5 x 10 4 K, V = 7500K, E { = -10 4 K, 
n el = 1.9, and A = 30K. 

FIG. 2. Temperature dependences of the linear susceptibility Xo- The pa- 
rameters are the same as in Fig. 1. The figure (b) shows the structures in 
low temperatures in (a). 

FIG. 3. Temperature dependences of the elastic constant c/c for g = 6.9K. 
The parameters are the same as in Fig. 1. In (a), the calculation is compared 
with the (en — Ci 2 )/2 constant of CeTe. In (b), the comparison with the c 44 
constant is made. Experimental data are shown by the dots. 
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